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Abstract. We derive a new formulation of the fully general relativistic equations describing a stationary equatorial 
MHD outflow from a rotating central object. The wind solution appears as a level contour of a "Bernoulli" function 
fixed by the requirements that it must pass through the slow and fa st critical points. This approach is the general 



relativistic extension to the classical treatment of Sakurai (1985). We discuss in details how the efficiency of 



the magnetic to kinetic energy conversion depends mainly on the geometry of the flux tubes and show that the 
magnetic acceleration can work very well under some conditions. We show how this tool can be used for the study 
of several astrophysical phenomena, among which gamma-ray bursts. 
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1. Introduction 

Magnetized winds are believed to be present in many as- 
trophysical objects. They were first put forward in the 
context of the solar wind. More recently they were dis- 
covered to probably play a major role in many situations 
where a relativistic flow is powered by a central rapidly 
rotating compact object: wind of pulsars, jets in radio 
galaxies, quasars, Seyfert galaxies and BL Lac objects, 
microquasars and even possibly gamma-ray bursts. 

The first quantitative model of a magnetic stellar wind 



the classical case. Goldreich & Julian (1970) studied cool 
isothermal relativistic winds. Kennel et al. (1983) ex- 



was developed by Weber & Davis (1967). They were con- 
sidering the equations of a stationary, axisymmetric, poly- 
tropic flow near the equatorial plane in classical MHD. 
They found that such a wind can carry off most of the 
angular momentum of the star and are very efficient to ac- 
celerate particles up to very high velocities. An important 
feature of magnetic winds is the existence of three "critical 
points" where the velocity of the flow equals the wave ve- 
locity of the three MHD wave modes (the slow, Alfven and 
fast modes) whereas in comparison non-magnetic winds 
(Parker 1958) have only one "critical point" where the 
velocity of the flow equals the sound speed (sonic waves 
being the only present wave mode). 

The first extension of this theory to relativistic winds 



is due to Michel (1969) in the context of radio pulsars. 
This work was considering cold outflows driven by rapidly 
rotating highly magnetized neutron stars. The main con- 
clusion was that the efficiency of the magnetic to ki- 
netic energy conversion was extremely low compared to 
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tended Michel's model to finite temperatures and relativis- 
tic injection speeds. All these works were always limited 
to the equatorial plane and either completely neglected 
the effect o f gravity or adop ted an approximative treat- 
ment for it. Okamoto (1978 ) first included an exact gen- 
eral relativistic description of the gravity fleld. His work 
was also not limited to the equatorial plane, applying for 
that the powerful concept of flux tubes. However he re- 
stricted h is study to pressur el ess flo ws only. In a series 
of papers Camenzind (1986a , b|, 1987) derived a complete 
set of equations describing a stationary axisymmetric rel- 
ativistic magnetic wind in an arbitrary metric. He then 
solved these equations in some particular cases (cold flows, 
jet geometries). 

The goal of this paper is to present a formulation of the 
equations governing a stationary axisymmetric MHD flow 
in the equatorial plane including an exact treatment of 
all effects (thermal pressure, gravity and arbitrary shapes 
of flux tubes) which allows a direct comparison with the 
classical model of Weber fc Davis (1967 ) , so that the rel- 
ativistic effects can be easily identifled. This is done in 
Sects. ^ and ||, where we worked by analogy with the for- 
mulation of the classical case by Sakurai (1985). Then we 
study in details the efficiency of the magnetic to energy 
conversion (Sect. in particular the influence of the flux 
tubes geometry and of the gravity. We conflrm and ex- 
tend the results of Begelman & Li (1994) and we show 



that a large variety of situations is expected from very in- 



efficient winds like those considered by Michel (1969) to 
highly efflcient cases. Because our model assumes axisym- 
metry and focus on the equatorial plane, it fully applies 
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only to simple astrophysical objects like isolated neutron 
stars. On the other hand it can also describe the outer 
parts of more complex objects, e.g. compact objects with 
accretion disks or complex magnetospheres, as long as the 
magnetic field can be approximated as monopole like at 
these distances from the source. In the particular case of 
gamma-ray bursts the possibility of Poynting-flux domi- 
nated fireballs is briefiy discussed in Sect. ||. This is sum- 
marized in Sect. ^. 

2. The wind equations 

The conservation laws of the general-relativistic magneto- 
hydrodynamics have been derived by Bekenstein & Oron 
(1978| ). pamenzind (1986t ) used these results to obtain 



the equations governing a stationary axisymmetric wind. 
In this section, wc first recall these equations and we ap- 
ply them to the particular case of a flow occurring in the 
equatorial plane. Then we derive a new formulation for 
this problem where the wind solution is a level contour of 
a Bernoulli-like function. A very similar formulation was 
earlier studied by Bakurai (1985 ) for the non-relativistic 
case. The similarity allows us to compare our results with 
those of the classical case. 



6^ is the magnetic field according to a comoving observer, 
written as 



6^ = F*"""! 



(6) 



with — bf^jD^^. In the comoving frame b^ reduces to the 
common magnetic field 

^'^^o = ^*^" = (0,S) . (7) 

All dissipative effects (heat conduction, viscosity, cooling 
by radiation, etc.) have been neglected so that the flow 
is adiabatic. Ideal MHD is also assumed, which means 
that the proper electric field as seen in the plasma frame 
vanishes 



F^'u^ = 



(8) 



The following quantities appear in these equations: u'^ is 
the 4-velocity, p is the comoving mass density, P is the 
pressure and h is the specific enthalpy, which is given by 
(assuming a constant adiabatic index 7) 



7 P 
7- 1 p 



(9) 



The last assumption is that the particle number (or mass) 
is conserved 



2.1. Assumptions and basic equations 

Any stationary and axisymmetric space-time can be rep- 
resented by the following metric: 



ds^ = ffttdt^ + 2gtcjidtd4) 



+ .gabdxMa;^ 



(1) 



where the coordinates t and cf) correspond to the two sym- 
metries of the space-time defined by the two Killing fields 
k — dt (stationarity) and m — (axisymmetry). The 
metric coefficients gab depend only on the two remaining 
coordinates (a — 1,2). Note that in the whole paper 
we will make use of the ( — hH — h) signature for the metric. 
The electromagnetic field is described by the field tensor 
Ff''' and the dual field tensor F*'"^ = ^e'^'P'' Fp„ {e^^P" 
being the Levi-Civita alternating tensor) which satisfy the 
Maxwell equations 



V^F*'"' = 



(2) 



The motion of the plasma is governed by the energy- and 
momentum conservation equation 



V^T^" = 



(3) 



where the energy- momentum tensor TP'^ is made up of the 
fiuid part 



and of the electromagnetic part 
1 52 



(4) 



(5) 



{puP) - 



(10) 



Before writing the wind equations, it is useful to define 
the specific angular momentum of the fiow 



(11) 



As a consequence of the symmetries, the non vanishing 
elements of the electromagnetic tensor are given by 



Fti 

Ft2 

F(f,i 

F^2 

F12 



-Fit 

-F2t 



cn{p) 



p 



-cn{p) 



viP) 
p 



-gu 



viP) 



~gu 



-Fi^ ^ -c 

'FlA, -C 
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viP) 
p 

viP)' 



-gu 



-gu 



F21 



P 



viP) 



-g{n{P)i 



(12) 

(13) 
(14) 
(15) 
(16) 



where the angular frequency of the streamline r2(P) at 
its footpoint and the mass flux per unit flux tube ri{P) 
are constant along each flow line P. The corresponding 
magnetic field is 



utu"- 



b' 
52 



viP) 

P 

viP) 

P Ut 
T](P) C2 
P Ut 

riiP) c2 



n{p) 
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n{p)i\ 2 



(17) 
(18) 
(19) 
(20) 
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It is also useful to give the expression of the classical mag- 
netic field in the frame of the central object for the com- 
parison with the classical case. In this frame, it is related 
to the dual field tensor by i?^ 



F,^o- Then 



B' = — 



1 



1 P 



^) gtt , 



B' = - 



B^ = 



{gtt + ^{P)gt^) V? 



c2 r](P) 
1 P 

c2 ^{p) 



(21) 
(22) 
(23) 
(24) 



Because the flovif is stationary and axisymmetric, the 
total angular momentum i(P) and the total energy 
Etot{P) are also conserved along each flow line P which 
provides us with two new equations 



L(P) = 
and 

EtotiP) 



Ut_ 
n2 



Amf{P) 



gtt + ^{P)gt4, 



I + (5t0 + n{P)g^4,) 



^^■\h + n{P) 



4717^2 (P) 



gtt + n{P)g 



't<j, 



I + {gt^ + n{P)g^^) 



(26) 



It is convenient to write the total energy as Etot{P) = 
+ E{P) + n{P)L{P) so that the energy conservation 



can have the simpler form: 



ut 



E{P) + = -h^ 



1 - 



n{p)i 



(27) 



Each flow line P is completely determined by the four 
constants 17(P), r]{P), L{P) and E{P). The hght surface 
is defined by 



gtt + 2gt^^{P)+g^^n\P)^Q 

and the Alfvcn point is fixed by two conditions 

4 [gtt + 2gt4,n{P) + g^^^\P)\A = ' 



1 {gt4> + ^{P)g<i><i>)t 



L{P) 



c2 Ml E{P) + c2 

where the "Mach" number M is given by 
ATirf{P)h 



M 



(28) 

(29) 
(30) 

(31) 



One sees immediately from ( p^ and (|2^) that the Alfvcn 
point stays always inside the light surface (because of 
Ma > 0). 



f. The 



2.2. The wind equations in the equatorial plane 

We use now the spherical coordinates (a;^ = r and x'^ 
and limit our study to the equatorial plane 6 
specification of the flow line P is dropped from here on 
and fi, rj, E, and L are used instead of ^{P), ?7(P), E{P), 
and L(P). Because of the symmetry, and B^ vanish 
in this plane but this is not necessarily the case for their 
derivatives. Then the conservation of mass ( |l0|) can be 
written 



gs{r)pu'^ ~ rh 



(32) 



where the function s{r) depends on the geometry of the 
flux tubes. In the simple case where dgu^ — and 
deB^ — 0, we have s{r) — const, (constant opening angle). 
Otherwise we have 



"^-■2 s(r) 

The conservation of angular momentum 
servation of energy (k^) read 



(33) 
51) and the con- 



(25) L 



h ^ 

gt^-^ + ^- 



- gttg4><t> $ p 

c2 A^niin? 

ft^- gttg.l><t, $2p 
c2 47rm2 



and 

E + 



= -h 



gtt + ^gtcji t gtrp + ^300 



(34) 



(35) 



where instead of using rj we have introduced the magnetic 
flux $ = rh/r]. The Eqs. (|2|), (||) and (|3|) are completed 
by the normalization of the four velocity 

gtt {u'Y + 2gt^u'u^ + g^^ (u^)' + g„ {u'^ f - -c^ (36) 

and the equation of state. We assume here, like in Sakuraj 
1(1985 ), a polytropic relation P = np'^ so that the specific 
enthalpy is given by 



/l = c2 + 



7- 



1 



up 



7-1 



(37) 



The system of Eqs. (jSg), (gj, (|35|), (g6|) and ^ de- 
scribes entirely the flow determined by the six constants 
f2, E^ L, $, m, K and the free function s{r). In addition the 
two supplementary conditions 

1 , 



{gtt 



2ngt^ + n'g^.^,)^ = Ml 



$2 /9AC2 



and 



.9*0 + 



gtt + 2flgt4, + fl'^c 



E 



(38) 



(39) 



(all quantities with index A are computed at the Alfven 
point) must be fulfilled, so that the flow remains regular 
at the Alfven point. 



4 



F. Daigne and G. Drenkhahn; Stationary equatorial MHD flows in general relativity 



The classical limit (for a weak gravitational field and 
for velocities small compared to the speed of light) of this 
system of equations in the case where s(r) = const. — 1 
gives exactly the Eqs. (1) to (6) in Sakurai (1985| ). Notice 
that E, $ and r2 have the same meaning in both papers 
whereas we use here different notations for the mass flux 
rh, the total angular momentum L and the polytropic con- 



stant K which are respectively /, llr^ and K in pakurai 



(1985|)] Notice also that the relation (|39D between L and 
TA tends towards L 
all notations are fully consistent 



fir^ in the classical limit, so that 



2.3. A Bernoulli-like formulation 



Following Bakurai (1985) we use the dimensionless vari- 
ables X — r/vA and y — p/pA- The metric coefficients 
are also normalized to become dimensionless gtt{x) — 
9tt/c^, 9t^{x) = gt^/{cr), g^^{x) = g^^/r'^, grr{x) = 



and 



-g/r^. These coefficients may be 



not only functions of x but also of some parameters 
defining the metric (they are given for the Minkowski, 
Schwarzschild and Kerr metric in appendix ^ . We define 

w'^ = {g^^ - gug4>4^ !(?■ We normalize s by s{x) = s/sa 

and define a dimensionless specific enthalpy hy h = h/c^. 
Concerning four-vectors like u^^ or we will use the def- 
initions i* = A*, i"^ = rsin6'yl'^, i'' = A'' and = rA^ 
so that the spatial part of these vectors is now given in the 
usual basis dr,^dg, j^^d^, where all components have 
the same dimension. We introduce four normalized pa- 
rameters 



e' 



1 



c2 y.'iAr\pA 
Ii^Pa 



{^taY 



(40) 

(41) 
(42) 
(43) 



and we are now able to rewrite the system of Eqs. (^ 
and dH) to (|7|): 



- - 2 ' 

g s X y- 



(44) 



c 



~ r Ml). 2 
- r atMi)~2 

gtcj>h + Vuj'j^zu xy 



A 



E' + 1 = -h 



gtt + vuj'gt4,x ] u 
u 



c 

2 I 

-1 = fftt ("*) + 2gt4,u* — + g^^ I — 



(45) 



(46) 



Hy) 



e 



,7-1 



r 



(47) 
(48) 



In ( ^6| ) and ( p5| ) the constants L/tac and <i>^pA/47rm^ 
have been eliminated using the conditions ( ^Sf ) and (39) 
at the Alfven point, which now read 



- (gtt(l) + 2V^5t^(l)+c^'g^^(l)) - Ml 



Attw? h{l) 

$2 PA 



and 

_ gt<^(l) -h\/u7g00(l) 
Ml 

so that 

£Va _ Ml) 
Ml 



L 



1 



TAC + 1 



47rm2 



and 



L _Ka , 
— -J^iE 



VAC 



1) 



(49) 
(50) 

(51) 
(52) 



with Ka = 5*0(1) + Vt^'g0</-(1)- 

The component 'ff' can be expressed from (^) and 
substituted into ( ^ ) to provide a first relation between w* 
and u'^. A second relation between these two components 
IS given by (^) when subtracting Ka/MIx (||). It al- 
lows us to express all components of the four velocity as 



functions of only x and y. Then the remaining Eq. (46) 
becomes a Bernoulli-like equation H{x, y) = const, as (10) 
of Sakurai (1985). We do not further elaborate on the dif- 
ferent steps which lead to this final expressions: 



Ky) = 1 + 



e' 



/IF? — A ^i^^y) 
V^(a;,y) 

D{x,y) 



- = VKix^y) 



c 



-gsx^y 



H{x,y) + l = vj'^xh{y) y/K{x,y) 
= E' + 1 , 



\Af{x,y)\ 



(53) 
(54) 

(55) 

(56) 

(57) 



where we have introduced the following auxiliary func- 
tions: 



K{x,y) 
N{x,y) 

D{x,y) 



= 1 



/3' 



-g s^x^y"^ 



Ma + V uj'Ka )gci>,t>x + KAgtct> Hy) 



-h{l)w'^xy 



(58) 



(59) 



^Ma + Vu?Ka ) gt<p x + KAgtt 
uj'H^)'^^x^y 



Hy) 
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N{x, y) = {jjtt + 2^/uJ'x gt^ + lo'x^ g^^^ H^)y 

+Mlh{y) , (60) 
1^{x,y) = - [g^^D'^{x,y) + 2gt4,D{x,y)N{x,y) 

+ ~guN^x,y)] . (61) 

Equation ( p^ ) is the Bernoulli equation we will now con- 
sider. The solution y{x) of the wind equations appears as 
the level contour E' of the surface H{x,y). Notice that 
the classical limits of O' and uj' are not exa ctly the 



corresponding /?, Q and uj parameters used by Bakurai 
(1985|]J who made the choice of normalizing these quan- 
tities with the value of the gravitational potential at the 
Alfven point G M/va [see Eqs. (11a), (lib) and (He) of 
Sakurai (198?: )] whereas we used to make the defini- 
tion of these parameters more general. However a simple 
relation applies between Sakurai's and our parameters: 



p/13' = e/e' = uj/uj' = va/i 



■g where Vg 



GM/c' 



IS 



the gravitational radius. This is also valid for the classical 
limit of the Bernoulli function H(x, y) and the definition 
used by Sakurai (1985 ). 

Before studying H{x, y) in the following section, we 
have to note that this function is not defined everywhere 
in the region x > 0, y > as it is the case in the clas- 
sical limit. The function 2?(x,?/) must be strictly positive 
so that H{x,y) is well defined and the velocities are not 
imaginary. The domain where this condition applies is de- 
termined in appendix^. In the sub-Alfvenic region {y > I) 
this domain lies always inside the light surface, its location 
is given more precisely in appendix 

3. Description of the solutions 

3.1. Properties of the Bernoulli function 

As described by Sakurai (1985| ) in the classical case, the 
Bernoulli function H{x, y) has the following properties: 

1. At ?; = 1 (p = Pa) H(x,y) diverges if x ^ I and 
remains finite if a; = 1 (r = va). It means that all 
solutions going from the sub-Alfvenic region (y > 1) 
to the super-Alfvenic region [y < 1) must pass through 
the Alfven point x = y = 1 which is the only "hole" in 
the infinite "wall" y = I. 

2. Two important curves in the x-y plane are the so called 
slow/fast mode Mach curve defined by 



dH 
dy 



{x,y)=0 



(62) 



(the slow mode corresponds to the sub-Alfvenic region 
y > 1 and the fast mode to the super-Alfvenic region 
y < 1) and the gravitational throat curve (so called by 
analogy with the de Laval nozzle) defined by 



(63) 



At the intersections of these two curves, the function 
H{x, y) is locally flat, corresponding to an X-type crit- 
ical point (or 0-type point). 




X = r/r^ 

Fig. 1. Solution plane of the v^rind equations. The 

gray region corresponds to the domain where the Bernoulli 
function is well defined. In the sub-Alfvenic region (y > 1), 
it is limited by the light surface. The thick line (Alfven 
mode) separates the sub- and the super-Alfvenic modes. 
The dashed line indicates the slow {y > 1) and fast {y < 1) 
mode Mach curves and the dotted line the gravitational 
throat curve. The gravitational throat curve is very close 
to the fast mode Mach curve for this particular case. The 
slow (S) and fast (F) critical points are the intersections of 
the Mach and throat curves. The Alfven point is indicated 
by A. 

3. All level contours of H{x, y) going from y — > -|-oo in the 
sub-Alfvenic region to y — > 0+ in the super-Alfvenic 
region must cross these critical lines. They have to 
cross them simultaneously to be not interrupted. This 
means that the solution must pass through two criti- 
cal points defined as the slow (respectively fast) crit- 
ical point Xs,j/s (resp. Xf,yf), intersection of the slow 
(resp. fast) mode Mach curve and the gravitational 
throat curve. This imposes two new conditions for the 
solution: 



Hixs,ys)=E' and H {x^yf) ^ E' 



(64) 



Figure ^ shows the x~y plane for a particular choice of 
the parameters with the slow/fast mode Mach curve, the 
gravitational throat curve. Different level contours of the 
function H{x, y) are shown in Fig. |^. The solution is one 
of the level contours and passes through the slow and fast 
critical points. These figures are very similar to Fig. 1 of 
Sakurai (1985|) obtain ed in the classical case. Notice that 
Bcgelman & Li (1994) restrict their study to cold flows. A 



consequence of this assumption is that a purely radial flow 
is a singular case where the fast critical point is at inflnite 
radius. The fast point moves inward to finite radii only if 
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Fig. 2. A few level contours of the Bernoulli function are 
shown. The physical solution (thick line) starts in the sub- 
Alfvenic region, crosses the slow mode Mach curve at the 
slow critical point (S) then reaches the Alfvcn point (A) 
and enters the super- Alfvenic region where it crosses the 
fast mode Mach curve at the fast critical point (F). This 
calculation has been made for a Schwarzschild black hole 
with ni — 0.01 and the parameters (7 = 4/3, Q' — 0.04 
and Lu' — 0.5) have been chosen so that the different points 
are well separated. 

the flow diverges over-radially, like in the situations we will 
study in the next section. This is not a necessary condition 
in our more general model. Because we include gravity and 
thermal pressure the fast point is always located at finite 
distances even in a purely radial flow. 



Fig. 3. Parameter space: there are no wind solutions in 
the lower left part ( "Static" : the limits are only approx- 
imatively indicated). The upper left part ("Centrifugal') 
corresponds to winds where the thermal pressure is negli- 
gible whereas in the lower right part ( "Thermal" ) the cen- 
trifugal force plays no essential role. For a Schwarzschild 
metric with m = 0.01 (so w' < 0.98 and oj < 98) we 
have computed several series of solutions with constant 
terminal Lorentz factor Too — 1.01, 1.5, 10, 100 and 1000 
(solid lines). For Too = 10 we show also the case where 
TO = 0.1 (dashed line: oj < 8) and to = 0.001 (dotted line: 
u! < 998). The dotted horizontal line corresponds to ^ = 1 
(initially the electromagnetic and the matter energy fluxes 
are equal) for m — 0.01. For to — 0.01 and Poo — 1-01 and 
1.5 we indicate at different positions the value of the effi- 
ciency ejf of the magnetic to kinetic energy conversion. 



3.2. Classification of wind solutions by dimension less 
parameters 



Equation ( |57| ) can formally be written as 



(65) 



where to stands for the parameters defining the metric 
(see appendix ^ for the definition of to in usual cases) 
and g stands for the parameters fixing the geometry of 
the flux tubes (i.e. deflning s{x)). We will for the moment 
restrict our study to the case where s — const. — 1 and 
we fix the adiabatic index to 7 = 4/3 for the rest of this 
paper. From the two conditions at the slow/fast critical 
points (|6^) we can fi x two parameters. In practice, fol- 
lowing Sakurai (1985), we adjust so that the value of 
H is the same at the two critical points and this value 
gives £". Therefore all solutions are determined by only 
three independent parameters: O', to' and to. Compared 



to the classical case studied by Sakurai (1985), there is one 



supplementary parameter: to. This is due to the presence 
of a characteristic length scale related to the structure of 
the space-time (typically the gravitational radius GM/c^) 
which has no classical counterpart. 

The signification of these three parameters is clear: to 
measures the intensity of the gravitational field, 8' gives 
the strength of the thermal pressure (the pressureless case 
which is often considered corresponds to O' = 0) and uj' 
measures the strength of the centrifugal force in acceler- 
ating the wind, or equivalently the effect of the magnetic 
field. 

In the limit of the flat space-time (Minkowski metric) 
it is well known that the slow point does not exist anymore 
and that the solution cannot start at an arbitrary small 
radius because it cannot cross the slow mode Mach curve 
without having an infinite derivative . We therefore cannot 
describe the region near the source but this is clearly be- 
cause gravity influences the solution signiflcantly at small 
radii and the Minkowski approximation breaks down. In 
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the case without gravity, as in the classical case, a sup- 
plementary parameter {/3') must be fixed (for instance by 
fixing the mass flux). In many cases the physical condi- 
tions at the basis of the wind are a complex and not well 
understood question, and some assumptions made every- 
where else in the flow are probably not valid here (like the 
adiabaticity) . In this context, just adopting a Minkowski 
metric and including in the /?' parameter all the unknown 
physics fixing the mass flux at the basis of the wind is more 
elegant than adopting a Schwarzschild (or Kerr) metric 
and applying the solution up to the source. 

Before exploring the three-parameters space we have 
defined, it is useful to express the relation between m, uj' 
and & and more usual physical quantities. If m is the ratio 
of the gravitational radius of the source over the Alfven 
radius, we have 



M 



(66) 



The angular frequency Q can be identified with the rota- 
tion rate of the source and fixes the value of lo': 



uj' = 0.11 • 



104 Hz 



106] 



(67) 



The local sound speed is defined by {(f/cf = "fP/ph so 
that 6' can be related to the ratio of the sound speed at 
the Alfven point (fp^ over the speed of hght (for 7 = 4/3): 



e' = 1.0 -10-2. f ^ 
0.1c 



i-3(cyc)^ 

0.97 



(68) 



The three parameters 0', lo' and m being fixed, the solu- 
tion of the wind equations is simply found by adjusting j3' 
so that the condition (^^ is fulfilled. For each step with /?' 
fixed, the slow and fast critical points are determined by 
a simple Newton-Raphson procedure. The exact expres- 
sions of H and its derivatives are used. We have explored 
in detail the parameter space and the results are presented 
in Fig. ^. Notice that for a given m, lo' is limited to the 
interval 



(69) 



where 



is due to the condition that the Alfven point 



lies inside the light surface and ti^^j^ is non vanishing only 
in the case of the Kerr metric. In this metric there are no 
solutions without rotation {lo' = 0) because the matter is 
forced to rotate in the vicinity of the central source. The 
analytical expressions of lo'^^^ and w^^x ^''^^ given in ap- 
pendix Here we have considered a Schwarzschild metric 
with different values of m. We show the results in Q-lo 
coordinates, where Q — Q' /m and lo — lo' /m are the 



parameters used by Sakurai (1985). Like in the classical 
case, there are no wind solutions in the lower left part of 
the plane ("Static") because neither the centrifugal force 
(magnetic acceleration) nor the thermal pressure are suf- 
ficient to power the wind. The limits of this region are 
only approximatively indicated. For a pure thermal wind 



{lo — 0) and a Schwarzschild metric, the minimal value of 
O is given by 



©min = (7 - 1' 



1 



(70) 



tending towards 7 — 1 for m ^ 0, in agreement with the 
classical case. In the pressureless case (8 = 0) the mini- 
mal value of LO tends towards (3/2)^/ ^ for to — > 0, which 
also corresponds to the limit given by Sakurai (1985 ). The 
upper left part corresponds to winds where the centrifugal 
force dominates and the lower right part corresponds to 
pure thermal winds. In Fig. ^ we have plotted several so- 
lutions with constant terminal Lorentz factor Too = 1.01, 
1.5, 10, 100 and 1000 for to = 0.01 and in the particular 
case Too — 10 we have also plotted the same curves for 
TO = 0.1 and TO = 0.01 to show the effect of varying the 
gravitational field. 

4. Efficiency of the magnetic to kinetic energy 
transfer 

4.1. Expressions of the energy fluxes 

To study the efficiency of the winds computed in the pre- 
vious section, we need to express the different components 
of the energy flux along the flow: 



-En 



57 



E + h- 



gttu +9t,j>— ) h{y) mc , 



(71) 



47rm^ 



TO flL 



= v Lo'xvo 7r- I yuj'xu 

Mi 



c 



E'total — En 



Eo 



m{c^ + E + VLL) 



{ 



E' 

gu{l) + 3(0(1) 



Ml 



(72) 



(73) 



Notice that the matter part is made up of the rest-mass, 
kinetic and internal energy of the matter. We define two 
parameters: the initial baryonic load rj 



1 -^matter 

77 



TOC^ 



(74) 



and the ratio of the initial power injected in the electro- 
magnetic field over the initial power injected in the matter 



Eo,- 



E„ 



(75) 
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where xq is the radius where the wind starts. The value of 
^ depends only weakly on xq , and we have arbitrary chosen 
Xq = 6m (or tq = Qrg). One sees that r] will be fixed by 9' 
[via the initial value of h{y)] whereas ^ depends strongly 
on uj' . Along the flow, the internal energy is converted 
into kinetic energy which accelerates the wind. Therefore 
if there were no magnetic field, the Lorentz factor at infin- 
ity would be Too = l/??- However, the magnetic field can 
also contribute to the acceleration (when coupled with the 
rotation) and depending on the efficiency of the conver- 
sion of electromagnetic into kinetic energy, the terminal 
Lorentz factor can be larger than l/?7, with a maximum 
value (complete conversion) given by 



-pmax 



(76) 



In reality the conversion will never be complete and will 
be estimated by the following fraction 



E, 

eff = l-- 



cm _ 1 



E, 



4.2. Inefficient conversion 



(77) 



In the pressureless case {&' — 0) it is well known that the 
magnetic to kinetic energy transfer is very inefficient for 
high terminal Lorentz factors ( [Michel 196E ) 



Too rhc^ _ J_ 

-E'total ^cxj 



(78) 



The curves for constant terminal Lorentz factor in 
Fig. 1^ show clearly that for highly relativistic winds {Toe = 
10, 100 and 1000) the terminal Lorentz factor Too is inde- 
pendent of UJ (or equivalently of ^) which means that there 
can only be a tiny magnetic to kinetic energy conversion. 
When UJ is very close to the maximal allowed value and 
the outflow is Poynting flux dominated (corresponding to 
the case where the Alfven point is at the light surface ra- 
dius) this tendency is not valid anymore. In this region the 
terminal Lorentz factor Too depends strongly on uj and is 
almost independent of (or equivalently of rj). However, 
even in this case only a tiny fraction of the magnetic en- 
ergy is converted into kinetic energy. The converted energy 
amount is great compared to the initial energy in the mat- 
ter part and therefore leads to a greater increase in F and 
-^matter throughout the flow. The efficiency eff of the con- 
version is maximal in the pressureless case (0 = 0) but is 
still rather small. In this case the parameters 77 and ^ are 
given by ~ 1 and ^ ~ -Etotai/wc^ — 1 = Fj^ — 1 which 
corresponds to 



eff 



1 



1 



FL 



(79) 



These tendencies are still present in mildly relativistic 
winds, as can be seen on Fig. |^ for F^o = 1.5 where we 
have indicated the evolution of eff along the curve. It is 



only within the classical limit that the conversion becomes 
important. This is shown in Fig. ^ [The case Fqo — 1.01 
{voo — 0.14 c)]. Here the efficiency of conversion reaches 
^ 60 % in the pressureless case, which is in agreement 
with the classical study of 3akurai (1985| ). 



4.3. Efficient conversion 

All wind solutions showing (except in the classical limit) a 
very bad efficiency of the electromagnetic to kinetic energy 
conversion have been calculated for a particular geome- 
try corresponding to s = const. = 1 in our dimensionless 
units. This corresponds to magnetic flux tubes of constant 
opening angle. Under the assumption that the velocity is 
purely radial and constant at inflnity it is possible to pre- 
dict analytically the asymptotic behavior of the flow for 
any kind of geometry s(x): 



y 



1 sx"^ 



so that the asymptotic expressions of the energy fluxes are 

-S'mattor /on^ 

(81) 



E, 



Uj'y^h{l) 1 1 



From the last equation, one sees that the magnetic to ki- 
netic energy conversion depends strongly on s. At infinity 
s — > is unphysical because it would mean that the energy 
diverges. The case where the opening angle is constant at 
infinity corresponds to E^ra ~* const. > at infinity so 
that the conversion is not complete and the case where the 
opening angle diverges (s +00) gives Ecm ^ so that 
the conversion is complete. These results indicate that all 
models considered in the previous section are inefficient 
due to a particular choice of the geometry: s = const. This 
assumption is certainly correct at very large distance from 
the source but the opening angle may have variations at 
smaller radii. Equation ( pT| ) indicates that every region 
where the opening angle increases is a region of efficient 
magnetic to kinetic energy transfer. This is in agreement 



with the results of Begelman & Li (1994). 

To check that the geometry is really the key param- 
eter governing the efficiency of such winds we have com- 
puted some models using various laws for the evolution of 
the opening angle s{x). The results are shown in Fig. ^ 
and confirm the previous analysis. We have considered a 
Schwarzschild metric with m = 0.01 and a wind model 
characterized by uj' — 0.97 and 6' = 0.1 (so the energy 
flux is initially dominated by the electromagnetic energy 
flux) . We plot the different energy fluxes and the "Lorentz 
factoi£j' as well as the geometrical function s{x) we have 

^ For commodity, we use the expression "Lorentz factor" for 
M* even if, strictly speaking, this should only be used at large 



F. Daigne and G. Drenkhahn: Stationary equatorial MHD flows in general relativity 



9 



1.4 
1.2 



6 
4 B- 



H 1 1 — I — h 



H 1 1 — h- 1- 



10 



Fig. 5. The influence of the geometry, gravitation and 
thermal energy upon the electromagnetic energy conver- 
sion. The solid, dotted and dashed lines correspond to dif- 
ferent gravitational field strengths, m — 0.001, m — 0.01 
and m = 0.1. a) The case for uj' — 0.95 oj'^^^^, O' — 
0.01 which is a cold Poynting-flux dominated outflow, 
b) Thermal energy dominated fireball: u' — O.luj'^^^, 
8' = 10. The dotted line lies very close to the solid line 
and is not visible, c) Non-relativistic case where the rest 
mass dominates with uj' = O.luj'^^^, O' = 0.1. 



angle increases, Eq. (|8^) shows that this relation should 
still be valid in the most general case, independent of the 
gravitational field or of the initial amount of thermal en- 
ergy. To check the validity of this result we consider 9 dif- 
ferent llj',6' ,m combinations, illustrating all possible situ- 
ations and for each of them we compute the evolution of 
the efficiency when varying Soc- As Fig. |^ shows, the ex- 
act shape of the geometry is not important, so we adopt 
a particular choice where s rises from sq = 1 to Soo be- 
tween xi = 100 and X2 = 200. This region lies always 
in the super- Alfvenic region, which as discussed above is 
the condition for magnetic to kinetic energy conversion. 
Figure |^ shows the quantity 



(83) 



E 






oo '^tx) 


E 


So 




Xq 



plotted over Soo/sq for the 9 different cases. Notice that 
with our choice of geometry s{x{) — sq. One sees that 
gravity and pressure changes the simple picture a bit. In 
the cold cases (Fig. ||a and Fig. jsjc) the converted energy 
fraction decreases for a stronger gravitational field and 
high values of m. On the other hand the gravitational field 
increases the energy conversion by a small amount in the 
hot thermal dominated case as seen in Fig. |b. But (H) 
remains still valid within a factor of 2. For the cases of low 
gravity and low thermal energy (solid lines in Fig. |^a and 
Fig. ||c) the quantity (|^ approaches 1 as expected. 

We can therefore conclude that the flow geometry al- 
ways dominates the energy conversion and all other pa- 
rameters play an only minor role. 



used in each case. Figure ^a corresponds to the inefficient 
case s = const. = 1. Figure ^ corresponds to the case 
where s increases in a region located between xi = 10 
and X2 = 18: the magnetic to kinetic energy conversion 
is immediately better. The efficiency eff increases also in 
geometries with different shapes (Fig. ||i,f,g,h) and dif- 
ferent locations of the s > 1-region, provided that this 
region lies beyond the fast point as shown by p3egelma"n 
& Li (: 994). In this the only relevant quantity 

which governs eff. A s > 1-region within the fast point 
like in Fig. does not increase ejj and is similar to the 
purely radial case. 

Begeiman & Li (1994) showed that the electromagnetic 
energy flux decreases like 



(82) 



for a cold flow in Minkowski metric. If the asymptotic 
regime is already reached in the region where the opening 





oo siXf) 







distance of the source where the metric is very close to the 
Minkowski metric, i.e. for x ^ 1. The correct expression of 
the Lorentz factor should be corrected with the lapse function, 
r = ai* , where a = 1/ -\/— — > 1 for a; — > -|-oo. 



5. Application to gamma-ray bursts 

Since the discovery of their optical afterglows gamma- 
ray bursts (hereafter GRBs) have been known to be lo- 
cated at cosmological distance. More than ten redshifts 
have been measured from z = 0.43 (GRB 990712) to 
z = 4.5 (GRB 000131). The corresponding radiated en- 
ergy in the gamma-ray domain (20 - 20000 keV) ranges 
from 5- 10^1 erg (GRB 970228) to 2-105'' erg (GRB 990123) 
assuming isotropic emission. Most sources that have been 
proposed to explain such a huge release of energy in a 
few seconds involve a rapidly rotating compact stellar- 
mass core. Among them the two most popular are merg- 
ers of compact objects (neutron stars binary or neutron 
star - black hole systems) or collapses of very massive 



stars to a black hole (coUapsars) (Meszaros & Rees 1992; 


Narayan et al. 1992; 


Mochkovitch et al. 1993 




Woosley 


1993 




Paczyhski 1998 


1. In both cases, the resulting sys- 



tem is a stellar mass black hole surrounded by a thick 
torus made of stellar debris or of infalling stellar mate- 
rial partially supported by centrifugal forces. An other in- 



teresting proposition (Usov 1992; Kluzniak & Ruderman 
1998; Spruit 1999| ) associates GRBs with highly magne- 



tized millisecond pulsars. The location of the detected op- 
tical counterparts, well inside their host galaxy and pos- 
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Fig. 4. Effect of the geometry on the efficiency of the magnetic to kinetic energy conversion: wc consider 
a Schwarzschild metric with m = 0.01 and a wind solution characterized by Q' = 0.65 and w' = 0.74. This corresponds 
to an initial energy flux which is dominated by the electromagnetic part: ^ = 3.1 for all models so that initially 

75 % of the energy is magnetic. All solutions presented here have 77 = 2. 10~^ except for cases e {rj = 1.7 10~^) 
and h {r] = 1.8 10~^). The slow and fast critical points are located at Xg = 3.0 10~^ and Xf = 2.0 except for cases 
d {xf = 1.6) and h {xf = 1.6). On each figure dlogs/dlogo;, s{x) and the different components of the energy flux 
(matter/cm and total) arc presented as functions of the radius x. The "Lorcntz factor" ?i* is also shown (dotted 
line). Three vertical dotted lines show the location of the slow (s), the Alfven (A) and the fast point (f). Case a: 
s = const. = 1. For this particular choice of the geometry, the conversion is extremely inefficient {eff = 3.8 10~^) 
and the terminal Lorentz factor equals Fqc = 50 = I/77. Case b: s increases between xi = 10 and X2 = 18 reaching 
a maximal slope dlogs/dlogx = 1. The efficiency improves a lot: eff = 0.26 and T^o = 90. Case c: same as b but 
s increases between x\ = 1000 and X2 = 1800 {X2/X-1 is the same). It changes neither the efficiency nor the terminal 
Lorentz factor. Case d: same as b but ,s increases between xi = 1.5 and X2 = 2.7 (x^jxx is the same), i.e. before the 
position of the fast point in the reference solution a. Again the efficiency eff = 0.24 and the terminal Lorentz factor 
Too = 88 are almost unchanged. Notice that the fast critical point has moved to be almost at x\. Case e: same as b 
but § increases between xi = 0.1 and X2 = 0.18 {x2/xi is the same), i.e. before the Alfven point. The efficiency is again 
very low: eff = 1.2 10~^ and = 58 ~ l/r]. Case J: same as b but with a maximal slope of dlogs/dlogx = 4. The 
efficiency is better: eff = 0.69 and Foo = 157. Case g: same as b but the region where s increases is larger: x\ = 10 
and X2 = 100. Again the efficiency is better: eff = 0.68 and Foo = 156. Case h: we have considered a case where s 
increases from xi = 0.1 to X2 = 10^ with a maximal slope dlogs/dloga; = 0.4. Almost 90 % of the magnetic energy 
is converted into kinetic energy {eff = 0.88) so that Fqo = 206. 
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sibly associated with star-forming regions, seems to fa- 
vor the coUapsar scenario. However the other propositions 
cannot be ruled out, at least for short bursts, for which 
no optical counterpart has been detected yet. 

Whatever the source is, the released energy must ini- 
tially be injected in a wind which eventually becomes rela- 
tivistic. The existence of such a relativistic wind has been 
directly inferred from the observations of radio scintilla- 



needed by a particle starting from the basis of a flow line 
to reach the Alfven point: 



tion in GRB 970508 ( [Frail et al. 1997| ) and is also needed to 
avoid photon-photon annihilation. The absence of signa- 
ture of this last process in the BATSE spectra of GRBs im- 
plies very high Lorentz factor for the wind: F ~ 100-1000 



(Goodman 1986; Baring 1995). The second step consists 
in the conversion of a fraction of the wind kinetic energy 
into gamma-rays, probably via the formation of shocks 



within the wind it self (|Rees fc Meszaros 1994| ; Paigne fc 



Mochk pvitch 1998]) . Such internal shocks are expected if 
the wind is generated with a highly non uniform distri- 
bution of the Lorentz factor so that rapid layers catch up 
with slower ones. In the last step, the wind is decelerated 
when it interacts with the environment of the source and 
the resulting external shock is responsible for the after- 
glow observed in X-ray, optical and radio bands. 

The origin of the relativistic wind is the most complex 
of the three steps in this scenario. Several proposals have 
been made but only few calculations have been performed 
so that none appears to be fully conclusive. However it 
is suspected that large magnetic fields play an important 



role. In a previous paper (Spruit et al. 2001) we have con- 
sidered different possible geometries of magnetic fields in 
GRB outflows and we have proposed that in many cases, 
dissipation of magnetic energy by reconnection should oc- 
cur. The model we have presented in this paper allows us 
to investigate these questions in more details. In particu- 
lar we focus on the case where the outflow generated by 
the central engine is initially Poynting flux dominated (in 
the following, we assume that only 10% of the energy is 
initially injected in the matter). To be consistent with the 
observations showing that at the beginning of the after- 
glow emission, the matter flow is highly relativistic, we 
also impose that the terminal Lorentz factor has a large 
value (in the following, we will adopt Fqo — 100). This 
implies a reasonable efficiency of the magnetic to kinetic 
energy conversion. The goal of the study presented in this 
section is to illustrate that there are geometries allowing 
such an efficiency and to discuss the possibility of mag- 
netic reconnection in this scenario. 



Spruit et al. (2001) have shown that for typical GRB 
outflows the MHD approximation is valid to very large 
distance (^ 10""^^ cm) which is the main assumption of our 
calculations. The second main assumption - the station- 
arity of the flow - is of course less justified in the case of 
GRBs. However we can estimate the time scale to reach 
the stationary regime in our wind solutions as the time 



1 



-dx 



(84) 



(in the source frame). Let us estimate this time scale in a 
particular case. We consider a Poynting flux dominated 
wind (we adopt ^ — 9.0 so that only 10% of the en- 
ergy flux is initially injected in the matter) with a mod- 
erately low initial baryonic load (we take 77 = 1/50). We 
impose that the terminal Lorentz factor is Fqo — 100. 
If there were no magnetic to kinetic energy conversion, 
the Lorentz factor at infinity would only be l/ij = 50. 
In order to get a final Lorentz factor of 100, we need 
to assume that the geometry allows an efficiency ejf = 
{v^oa - 1) /C = 1/9- We have shown in Sect. § that this 
implies Soo/sq ~ ^/ (1 -f ^ — rjToo) = 1.125. For a given m, 
the two other parameters Q' and w' are fixed by the values 
of ^ and 77. We find that the following set of parameters: 
TO = 0.069, e' = 3.8 and uj' = 0.78 fulfiU the requirements 
and corresponds to a reasonable value of the Alfven radius 
TA and the angular frequency Q in the case of a millisec- 
ond pulsar-like source (A/ — I.4M0) which is most likely 
leading to an equatorial fiow as we are considering here: 
TA = 3.0 • 10^ cm and fl^S.S- 10^ Hz. Figure | shows the 
evolution of the "Lorentz factor" and the electromagnetic 
and matter energy fiuxes in this case. We have assumed a 
simple geometry like those in Sect. |^ where s increases in 
a region located between xi = 300 and X2 = 900, well out- 
side the fast critical point radius. The corresponding time 
scale to reach the stationary regime tstat is between ta/c 
and ~ 2rA/c, depending on the adopted value of the initial 
radius a;o. As ta/c = 10~^cm/s here, this is compatible 
with the timescale of the variability observed in GRBs 
profiles. This means that when the physical conditions at 
the basis of the flow vary on a time scale tvar ^ 1 ms 
the flow reacts instantaneously to reach a new stationary 
state corresponding to the new boundary conditions. Thus 
our calculation is a good approximation for the relativistic 
wind of GRBs. If the wind produced by the source lasts 
for a duration iw, our solution is appropriate for the phys- 
ical quantities within the corresponding shell when it is 
located at radius r. 

On the solution we present on Fig. |6[ the accelera- 
tion occurs in two phases. First the initial thermal energy 
is converted into kinetic energy, the magnetic energy re- 
maining unchanged. This phase ends at r ~ 10^ cm where 
Foo — 1/?? — 50. The second phase occurs in the region 
where the opening angle increases. Here a magnetic to ki- 
netic energy conversion takes place. We define the accel- 
eration radius race as the radius where the fiow reaches a 
Lorentz factor of F = 0.95F(x; and the acceleration can be 
considered as finished. The value of this radius is com- 
pletely dominated by the unknown fiow geometry and 
equals race — 2.7 • 10^ cm in this case. Even if the location 
of the region where the opening angle diverges would ex- 
tend to higher radii up to 10^°-10^^ cm, this radius is well 
below two other important radii: the photosphere radius 
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Tph where the wind becomes transparent and the recon- 
nection radius r^cc where the reconnection of the magnetic 
field should occur. These two radii have been estimated in 



Spruit ct al. (2001). The photosphere radius is the solution 
of 



2f2 



dr 



(85) 



and is independent of the duration of the burst as long 

(r/100)^(tw/ls). Here 

-,51, 



as rph < 2r2ciw = 6 • 10^'' cm 

we have rph = 6.2 • 10^° cm for i?totai = 10^^ erg/s and 
r-ph — 5.9 • 10^^ cm for ii^totai ~ 10^^ erg/s. This interval is 



marked by a thick line on Fig. || As the remaining thermal 
energy in the wind at such a large radius is very small, 
our adiabatic wind solution applies up to the reconnection 
radius, where magnetic dissipation starts. This radius is 
given by 



1/2 



(86) 



where e < 1 is a numerical factor of order unity measur- 
ing the reconnection speed in unit of the Alfven speed. In 
our case we have riec ~ 1.1 • 10^^ cm • e. As the magnetic 
energy flux is still 80% of the total energy flux at Tree, a 
very large amount of energy can possibly be dissipated at 
this large distance. Depending on the value of ii'totai and e, 
such reconnection events may start when the wind is still 
optically thick (low e, high E'totai) or when the wind is 
already transparent (high e, low -Etotai)- As the dissipated 
magnetic energy is probably first converted into thermal 
energy, the consequences for the wind may be very differ- 
ent in these two cases, (i) if the wind is optically thick, 
this injection of thermal energy should be converted, at 
least partially (up to the photosphere radius) into kinetic 
energy, leading to a third phase of acceleration; (ii) on the 
other hand, if the wind is transparent, reconnection events 
could directly contribute to the observed emission. Notice 
that all the radii we have computed are usually small com- 
pared to the typical radius where internal shocks occur 
(with Too = 100) 



ris ~ 3 • lOi 



Is 



(87) 



where ^var is the typical time scale of the variability in the 
initial distribution of the Lorentz factor and also small 
compared to the deceleration radius where the external 
shock becomes efficient (with Too = 100) 



rdc 



5 • 10^^ cm • 



E 



2052 erg 



1/3 



V 1 cm / 



-1/3 



where n is the density of the external medium and E the 
total energy of the wind at this radius. So these two "stan- 
dard" mechanisms are not affected by the reconnection 
events. However the relevant energy flux will be the kinetic 
energy flux at race, possibly increased to a larger value if 
the reconnection starts in the optically thick regime. 




Fig. 6. Geometry, "Lorentz factor" and energy fluxes for 
our example. The vertical dotted lines mark the radii of 
the slow-, the Alfven-, fast point and the acceleration ra- 
dius. 



6. Conclusions 

We have presented here a new formulation of the equa- 
tions governing a stationary axisymmetric MHD flow in 
the equatorial plane. This formulation includes an exact 
treatment of all effects: thermal pressure, gravity and ar- 
bitrary shape of the magnetic flux tubes. The wind solu- 
tion appears as the level contour of a Bernoulli-function 
which passes through two particular points: the slow and 
fast critical points. It allows a direct comparison with the 
classical model of Weber fc Davis (1967[ ), in particular in 
the formulation given by Sakurai (1985). Thus the specif- 
ically relativistic effects are easily identified. 

We have used our model to extend the study of the 



magnetic t o kinetic energy conversion made by Begelman 
|fc Li (1994|) . We show that the main parameter which fixes 
this efficiency is the shape of the magnetic fiux tubes. In 
the case of a constant opening angle, non-relativistic flows 
have a good efficiency of the magnetic to kinetic energy 
conversion but as soon as the terminal Lorentz factor is 
greater than ^1.5, this efficiency decreases rapidly. Such 
relativistic winds are not able to transfer a large fraction of 
their magnetic energy to the matter. On the other hand, 
regions where the opening angle diverges from the con- 
stant case are very efficient in converting magnetic into 
kinetic energy, even in the ultra-relativistic case. This is 
true as long as such regions are located beyond the fast 
critical point. Gravity and the thermal pressure play only 
a minor role. 

In Sect. 1^, we apply this model in the context of 
gamma-ray bursts (GRBs hereafter). In the case where the 
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wind produced by the source of GRBs is initially Poynting 
flux dominated, we have shown that the efficiency of the 
acceleration strongly depends on the geometry of the mag- 
netic flux tubes. We found that a large variety of situations 
is expected. If the magnetic tubes have the possibility to 
diverge strongly from a constant opening angle, it is pos- 
sible that most of the energy is eventually in kinetic form. 
On the other hand it is very likely that the magnetic to 
kinetic energy conversion is incomplete and that the wind 
is still Poynting flux dominated when it has reached its 
terminal Lorentz factor. We have demonstrated on one 
example that such a wind can lead to very promising sit- 
uations compared to the standard picture: a large amount 
of the magnetic energy can be dissipated at large radii by 
reconnection. This reconnection can start when the wind 
is optically thick or already transparent. So the large mag- 
netic energy reservoir could have two effects: a supplemen- 
tary acceleration phase increasing the final magnetic to 
kinetic energy conversion efficiency and/or a direct con- 
tribution to the emission. These two possibilities will be 
investigated in a future work. 
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Appendix A: Metric coefFicients 

The following table gives the metric coefficients in 
normalized units in three cases: the Minkowski (M), 
Schwarzschild (S) and Kerr (K) space-times. 





M 


s 


K 


gtt 

9t<i> 
9<t><t' 

9rr 


-1 



1 

1 

1 
1 


-1-1- ^ 

X 



1 

X 

1 


1 1 2m 

5 

1 _ 2m ^^2rn. 

X x^ 

1 






2m 


m (1 + Vl - a?-) 


Xq 






2m 










2am 
l+a^m^+2a^m-' 


\/ ^max 


1 


1 - 2m 


2am^ + \/l — 2m-t-a^m^ 
l+a2m2+2a2m^ 



where m = GMjr^f? and a — Jc/GAP (where J is the to- 
tal angular momentum of the black hole and < a < 1). 
The radii and Xc are respectively the radius of the 
horizon and of the ergosphere. We consider only the case 
where < m < ^ (the Alfvcn point is outside the ergo- 
sphere) . The minimum value w^jjj of lo' corresponds to the 
condition > (positive total angular momentum L) 
and the maximum value w^^x corresponds to the condi- 
tion M\ > (the Alfvcn point must be inside the light 
surface). 
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Appendix B: The light surface 

The light surface is defined by ( p8| ) which limits the region 
defined by 

(^{x) = a^m^uj' -1 + 2— (l- amVuyy +uj'x'^ 

< . (B.l) 

In the general case (with uj' < Wmaxi this corresponds to 
a domain < ij" < a; < x^ including the Alfven point 
{xa = 1). In the Minkowski case, x^^ — and x^^ = 

Appendix C: Domain of definition of the Bernoulli 
function 

The Bernoulli function H (x, y) is defined for I?(a;, y) > 
which gives the following condition: 

A{x)Y^ -2B{x)Y + C{x)>0 (C.l) 

where Y = h{\)y / h{y) is a function of y and O' only and 
A, B and C are functions of x, cj', a and m only (notice 
that it is completely independent of the function s{x)). 
The function Y{y) is strictly increasing from Y{Q) = to 
F(-l-oo) = -\-oo with Y{\) — 1. So we can focus to the 
pressureless case 0' = and Y — y, all other cases 8' > 
corresponding only to a contraction of the domain along 
the y-axis. The coefficients A, B and C are given by 

A{x) = -{■&j''xf ('{x) (C.2) 
B{x) = Ml{w^xf (C.3) 
C{x) = {w^{l)^^{x) (C.4) 
- A^A {9tti^)94,4,x^ + '2gt4,{^)gt,i,x + g4,4,{l)gtt) } 

The radius x being fixed, the equation ( |C.l| ) has 1 or 2 
positive roots delimiting the domain where (x, y) is well 
defined. Many configurations are possible for a Kerr metric 
and we do not specify them in details here. We discuss 
only the case of the Schwarzschild metric. Four different 
configurations are possible depending on m (with a critical 
case at m, = 1/3) and w' (with a critical case at lo' = 
27w? (1 - 2m)^). This is illustrated on Fig. |c[ 
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